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SYMBOLS 


H(s) Laplace transfer function of continuous filter 
H(z) Z-transform of H(s) 

j imaginary operator, 

n 0) /w , normalized sample frequency or discrete time index 

s n 

s Laplace variable 

T sample period, sec 

U(s) Laplace transform of u 

u continuous input variable 

Y(s) Laplace transform of y 

y continuous output variable 

Z Z-transform operator 

z Z-transform variable 

? damping ratio of second-order model 

T time constant , sec 

0) frequency of sinusoidal input variable, rad/sec 

0)^ corner frequency of first-order high-pass filter, rad/sec 

umdamped natural frequency of second-order model, rad/sec 

0)^ frequency of computation or sampling, rad/sec 

(") time derivative 

I function evaluated at s = 

' s = 
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SUMMARY 


The fact that aircraft motion simulators inherently have less travel than 
the actual aircraft being simulated has resulted in the application of digital 
washout filters in the computer software. Four commonly used mathematical 
models of linear first- and second-order high-pass washout filters were ana- 
lyzed. These models were Euler's Integration, Zero-Order Hold, Bilinear 
Transformation, and Second-Order Adams-Bashforth Integration. Bode responses 
for each model at various sample rates were compared to the continuous filter 
response. The results show that higher sample rates produce Bode responses 
approaching the continuous response and the Bilinear Transformation model pro- 
duced the best responses over the frequency spectrum and sample rates. Pole 
location analysis of each model in the z-plane shows the Bilinear Transforma- 
tion and Zero-Order Hold models gave stable poles regardless of time step 
size, whereas the other models did not always display stable poles. 

A near constant gain error over the entire frequency spectrum was dis- 
covered in the Zero-Order Hold cases and a correction gain was calculated for 
the first-order high-pass filter case. 


INTRODUCTION 


One of the continuing areas of research in real-time man- in- the- loop 
aircraft simulations is simulator motion recovery. The fact that motion simu- 
lators have less travel than the aircraft being simulated results in compro- 
mises in the motion felt by the simulator pilot. In an attempt to utilize the 
limited travel most effectively, low-frequency portions of the motion signals 
are eliminated with "washout filters." 

Washout filters are merely high-pass filters and can be implemented 
digitally by any number of mathematical models which approximate continuous 
filters. This study examines the effects of modeling a linear first- and 
second-order high-pass filter with four commonly used models. The models are: 
Euler's Integration, Zero-Order Hold (ZOH) , Bilinear Transformation (Tustin's 
Approximation) , and Second-Order Adams-Bashforth Integration. 

As a means of comparison, the frequency responses for each mathematical 
model are analyzed against the continuous filter response while varying the 
sample rate. The z-plane pole locations of each model are also analyzed for 
stability. 


HIGH-PASS FILTER MODELS 


To obtain the Bode plots for comparing the various mathematical models, 
the mathematical representations of the first- and second-order high-pass 
filters are needed. 


First-Order High-Pass Filter 

The first-order high-pass filter (FOHP) is defined as: 

• , 1 _ • 
y + - y = u 


( 1 ) 


Taking the Laplace transform of equation (1) and assuming all initial condi- 
tions are zero, one obtains the first-order transfer function: 


Y(s) 

U(s) 


H(s) 


s 


S + 0) 

c 


where u) 

c 


1/t is the cutoff frequency. 


( 2 ) 


Second-Order High-Pass Filter 
Defining the second-order high-pass filter (SOHP) as: 

y + 2^o)^y + (jj^y = u (3) 

and taking the Laplace transform of equation (3) (assuming all initial condi- 
tions are zero) yield the second-order transfer function 

Y(s) - 

U (s) + 2 ^ 0 ) s + 

^ n n 


MATHEMATICAL MODELS 


The mathematical models to be considered are Euler's Integration, Zero- 
Order Hold, Bilinear Transformation, and Second-Order Adams-Bashforth 
Integration. 


Euler's Integration 
Euler's Integration is defined as: 

y(n + 1) = y(n) + Ty(n) 


(5) 


2 



First-order filter with Euler implementation- Applying equation (5) to 
(1) and taking the z-transform of the equation result in the transfer 
function: 


H(z) 


z - 1 

z + u T - 
n 


1 


( 6 ) 


Second-order filter with Euler implementation- The transfer function for 
the second-order high-pass filter (equation (3)) can be written as two differ- 
ential equations: 


X = -2CW X - 03^X + u 

1 n 1 n 2 


X = X 
2 1 


(7) 


with the output equation as: 


y = -2?to X - (jj^x + u 
n 1 n 2 


Applying Euler's Integration and taking the z-transform, one obtains: 


H(z) = 


Z^ + Aj^Z + Ag 
z^ + Bxz + Bq 


( 8 ) 


(9) 


where 


Ao = 1 


A = -2 
1 


B = 1 - 2i;u) T + 

0 n n 

B^ = 2(?o)J - 1) 


Zero-Order Hold 

The z-transform of the ZOH implementation can be stated as: 


H(z) = ^ Z 




( 10 ) 


FOEP with ZOH implementation- Substituting equation (2) into (10), the 
z-transform of the FOHP becomes: 


H(z) 


z - 1 


z 


-u.„T 


( 11 ) 
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SOHP with ZOH implementation- Applying equation (10) to (4) yields: 


where 


H(z) 


+ Aj^z + A(j 

z^ + B^z + Bg 



\ = -a + Bo> 



Bj^ = -2e ^ cos(m) 

G = cos(k) + (m/k)sin(k) 


k = 40 ) T 
n 

i = (o) T)^ 
n 

m =0) t/1 - 42 
n 


( 12 ) 


Bilinear Transformation 

The z-transform for the Bilinear Transformation is defined as: 


H(z) 


H(s) I 
' s 


2 z - 1 

T z + 1 


(13) 


FOHP using bilinear- Applying equation (13) to (2), we obtain the 
z-transform of the FOHP as: 


H(z) 


2(z - 1) 

(to T + 2)z + (o) T - 2) 
c c 


(14) 


SOHP using bilinear- Applying equation (13) to (4), one finds the 
z-transform of the SOHP to be: 


AjZ^ + A, z + Aq 

H(z) = 

B,z^ + B^z + B^ 


(15) 


4 


where 



B. = 4 - 4^0) T + (o3 T)^ 

0 n n 

B, = -8 + 2(o) T)^ 

1 n 

B = 4 + 4?o) T + (o) T)^ 
^ n n 


Second-Order Adams-Bashforth Integration 
The Second-Order Adams-Bashforth Integration is defined as: 


y(n + 1) = y(n) + [3y(n) - y(n - 1)] 


(16) 


The z-transform of equation (16) is: 

H(z) = H(s) 1 


2z 


z - 1 
3z - 1 


(17) 


?0EP using Adams -Bash forth Integration- If we substitute equation (2) 
into (17), we find the FOHP transfer function as: 

2z(z - 1) 


H(z) = — 5- 

2z^ + (3o) T - 2)z - w T 


(18) 


SOHP using Adams-Bashforth Integration- Substitution of equation (4) 
into (17) results in the SOHP transfer function: 


H(z) = 


z^(A^z^ + A^z + Ajj) 

B4Z** + B3Z* + B2Z^ + B^z + B^ 


(19) 


where 



= (o) T)2 
0 ' n 

B, = 2[2Cco T - 3(o) T)^] 
n n 
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B = 4 - 16zui T + (3w T)^ 
n n 

B, = 4(3^0) T - 2) 

^ n 

B, = 4 

The z-transfer functions of the four mathematical models are summarized 
in table 1. Bode plots for each model were produced by evaluating these 
transfer functions at z = e^^"^ = cos(ojT) + j sin(mT). The amplitude and 
phase of the steady-state response to a sinusoidal input are plotted as a 
function of normalized frequency (normalized to the natural frequency of the 
filter) for various normalized sample frequencies (parameters with units of 
samples per cycle at the natural frequency of the filter) . The Bode responses 
of the associated continuous high-pass filters are included in the plots for 
comparison. 

The normalized sample frequency is defined as: 


for the FOHP, and 


for the SOHP. 

In the pole- location analyses, the poles of each model in table 1 are 
plotted as a function of o)j,T for the FOHP cases and w^^T for the SOHP 
cases. For the SOHP cases, damping ratios of 0.3, 0.707, and 0.9 are repre- 
sentative of purposes of these analyses. 


n = 


2tt 

Too 


n = 


n 


2tt 

To) 

n 


RESULTS 


First-Order High-Pass Filter 

Figures 1 through 4 present the Bode plots for the FOHP models. Based on 
the deviation from the desired continuous frequency response, the application 
of the Bilinear Transformation produced the best results over the sample fre- 
quency domain. All the mathematical models produced acceptable results when 
the normalized sample frequency (n) was greater than 50, or when the natural 
frequency times the sample period is less than or equal to 0.125 
[o}j,T < (2 tt/ 50)]. With the normalized sample frequency greater than 50, the 
gain deviation was within 1 dB and the phase was within 2° from the continuous 
FOHP for all mathematical models in the entire normalized frequency range up 
to the normalized Nyquist frequency of co/co^, = n/2. 
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Figures 5 through 8 show the pole locations of each model as a function 
of Euler's Integration (fig. 5) becomes unstable when cJ(,T > 2.0. 

The Second-Order Adams-Bashforth (fig. 8) becomes unstable when > 1.0, 

which corresponds to n < 6. This is also shown in the phase plot of fig- 
ure 4. In the limit, as w^,T approaches infinity, the Bilinear pole goes to 
-1.0 and the ZOH pole goes to 0, which means that these two models give stable 
poles over the range of u^T. Table 2 summarizes the value of 0)^,1 for mar- 
ginally stable poles. These results are also applicable to first-order low- 
pass filter models. 


The Bode plot for ZOH (fig. 2) shows a near-constant gain error across 
the frequency range for each sample rate. From the plot, one concludes that a 
DC gain must be calculated and applied to eliminate the error when using this 
model. The gain error is analyzed in appendix A and is found to be approxi- 
mately equal to: 


gain error = 



2rr 

n 



- 1/2 


Therefore, the ZOH model for FOHP should be: 


H(z) = 


z - 1 


(gain error) (z - e 


-0) T 
c 


Second-Order High-Pass Filter 

Comparing figures 9 through 20 with the addition of the damping ratio (?) 
as another variable to the Bode plots, the Bilinear Transformation again pro- 
duced the best overall results. With the normalized sample frequency greater 
than 50, the gain deviation was within 2.5 dB of the continuous SOHP for all 
mathematical models over the frequency range up to and including the appro- 
priate Nyqulst frequency. The phase deviation was within 9° of the continuous 
SOHP for all mathematical models except the ZOH model. However, the ZOH model 
was within 9 ° of the continuous SOHP when the normalized frequency was greater 
than 0.4. The ZOH results in figures 10, 14, and 18 show a gain adjustment 
can be calculated to compensate for the near-constant gain errors across the 
frequency spectrum, but this would not improve the phase errors. 

Figures 21 through 24 show the pole locations of all four models versus 
WnT. Euler's Integration (fig. 21) and Adams-Bashforth (fig. 24) show at what 
values of w^T the roots become unstable. These values of oj^T relate to 
the normalized sample frequency (n) and to the increased phase deviations in 
the associated Bode plots. Table 3 summarizes the values of co^T for mar- 
ginally stable poles. These results are also applicable to second-order low- 
pass filter models. Appendix B develops a stability check for Euler's 
Integration. Stable roots for SOHP using Euler's Integration will result if 

ojj^T < 2?. 
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CONCLUSION 


The application of all four mathematical models to the linear first- and 
second-order high-pass filters reflects the dependence of the digital models 
on sample rate. As seen from the Bode plots, reasonable results can be 
obtained when the normalized sample frequency (n) is greater than 50. This 
sample frequency produced gain deviations within 1 dB and phase deviations 
within 2° of the associated continuous first-order high-pass filter output to 
a sinusoidal input. Gain deviations within 2.5 dB and phase deviations within 
9° were obtained for the second-order high-pass filter except for the ZOH 
model at low normalized frequencies. The pole- location analyses show the 
Bilinear Transformation and ZOH models yield stable poles for all practical 
natural frequencies and calculation rates. Based on frequency response and 
pole-location analyses, Euler's Integration and Adams-Bashforth models should 
both be used with caution. The Bilinear Transformation model gave the best 
overall results over the entire frequency range up to the Nyquist frequency 
and is recommended for both high-pass filter cases. 

Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 94035, September 15, 1980 
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APPENDIX A 


APPARENT GAIN ON ZOH FIRST-ORDER HIGH-PASS FILTER 


The Bode plot for the FOHP using ZOH (fig. 2) shows a near constant gain 
across the frequency range as a function of sample rate. The following ana- 
lysis is an effort to find a compensation for the apparent gain error. 


ZOH Model 

The transfer function for the ZOH first-order high-pass filter model is: 

2-1 


H(z) = 


-0) T 
c 


z - e 

The frequency response equation is obtained by substituting the relationship; 




j/2iTa)\ , /2 ttu)\ 

= e-' I I = cos I I + 1 sin I I 


in the transfer function. If a = 2Tr/n and b = co/w , the equation is: 


juT. _ cos(ab) - 1 + j sin 


(ab) 


The magnitude squared is: 


H 


cos(ab) - e + j sin(ab) 


2[1 - cos(ab)] 

1 - 2e ^ cos(ab) + e 


Taking the series expansion of each term to the fifth power and combining the 
terms result in the equation: 


H 


b^[l - (bVl2)] 


(1 - a)(l + b^) + (7/12)a^ + a“b='{ (1/2) [ 1 - (h^/6)] - (a/6)[l - ih^/2)]} 
If b < 1 , the magnitude squared is : 


H 


(1 - a) + (7/12)a^ + (a‘‘b^/12)[6 - b^ - a(2 - b^) ] 
Therefore, for small b, the magnitude is 


. 1/2 


H = 


a - a + (7/12)a^ 
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Continuous Case 


The transfer function for the continuous case is: 


H(s) = 


TS + 1 


s = jo)T 


and the frequency response becomes 


H(jmT) = 


j (wT/m^) 
j(toT/o)^) + 1 


The magnitude is: 


jb 

jb + 1 


b" + 1 


/b^ + 1 


Therefore, for b « 1, 


H 1 b 


Comparing the magnitudes of the ZOH and continuous case results in: 


gam error 


n 12 \n / 


The gain error is a function of the normalized sample frequency (n) , if nor 
malized frequency (b) is less than unity. Therefore, to correct the gain 
error in the ZOH FOHP, the transfer function should be: 


H(z) = - 1 ) 

-w T 
c 

z - e 


gain = 


gam error 
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APPENDIX B 


EULER'S INTEGRATION - STABILITY CHECK FOR SOHP 


This appendix derives a simple check for stability when applying Euler's 
Integration on a SOHP. The transfer function for SOHP using Euler's Integra- 
tion is: 

_ - 2z + 1 

z^ + 2(?ai T - l)z + (1 - 2Ccj T + co^T^) 
n n n 

denominator is : 

- 1) + iu) Tv^l - ][z + (CO) T - 1) - jo) T/1 - 

jn -'n 

For stability, the magnitude of the vector Z must be less than or equal to 
unity (|zl < 1). Therefore, 

, , 1/2 

[(CO) T - 1)^ + (ui T/1 - ?2) ] < 1 

n n ~ 

which reduces to oi^^T < 2^ for stable roots. This relationship is also 
applicable when applying Euler's Integration to a second-order low-pass filter 
because the characteristic equation for both filters is the same. 


H(z) 

After factoring, the 

[z + (?a) T 
n 
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TABLE 1.- Z-TRANSFER FUNCTIONS OF THE FOUR MODELS 



d) 

u 

d) 


•H 

W 

/— N 


Vw' 

CO 

o 

o 

It 

o 


\ 


M 

♦O' 

I 



H 

3 

II 
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TABLE 2.- FOHP w T VALUE FOR MARGINALLY STABLE POLES 
c 


Math model 

H 

U 

3 

Euler 

ZOH 

Bilinear 

Adams-Bashforth 

2.0 

None 

Infinity 

1.0 


TABLE 3.- SOHP o) T VALUES FOR MARGINALLY STABLE POLES 
c 


Math model 

0) T 
n 

C = 0.3 



Euler‘S 

0.6 

1.41 

1.8 

ZOH 

None 

None 

None 

Bilinear 

Infinity 

Infinity 

Infinity 

Adams-Bashforth 

0.84 

0.95 

0.99 


“^Stable poles for Euler's Integration can be cal- 
culated using the relationship to^T < 2?. 


















Figure 5 



- Euler's integration FOHP pole location 



Figure 6.- Zero-order hold FOHP pole location. 




Figure 7.- Bilinear Transformation FOHP pole location. 



Figure 8.- Second-Order Adams-Bashforth FOHP pole location. 
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